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o : 

. We confront dark energy models which are currently similar to ACDM theory with observational 

O ' data which include the SNe data, matter density perturbations and baryon acoustic oscillations data. 

' DE cosmology under consideration may evolve to Big Rip, type II or type III future singularity, or to 

, Little Rip or Pseudo-Rip universe. It is shown that matter perturbations data define more precisely 

the possible deviation from ACDM model than consideration of SNe data only. The combined 
data analysis proves that DE models under consideration are as consistent as ACDM model. We 
demonstrate that growth of matter density perturbations may occur at sufficiently small background 
I I ■ density but still before the possible disintegration of bound objects (like clusters of gala^xies, galaxies, 

O ' etc) in Big Rip, type III singularity, Little Rip or Pseudo-Rip universe. This new effect may bring 

the future universe to chaotic state well before disintegration or Rip. 

^ ■ 

' ^! I. INTRODUCTION 

(N : 

' The late-time cosmic acceleration which seems to be proved by astrophysical data (see refs.[l|, 0]) opened the door 
for a number of (often exotic) theoretical models of so-called dark energy (DE) (for recent review, see 0, The 
dark energy contributing of nearly 72% of the total mass energy of the universe [5| has quite unusual properties like 
negative pressure and/or negative entropy, invisibility in the early universe, non-coupling with baryonic matter and 
effective non-observability, etc. The DE properties may vary in wide limits, what depends from the specific DE model 
under discussion. 

The negative pressure leads to that the equation of state (EoS) parameter ifDo is negative: 
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Wo = pdo/pdo < , (1) 



k>( ■ where pD is dark fluid energy-density and is the pressure. The subscript denotes that related quantities are 
j_j \ considered at present time. 

d • The data of observations (supernova data, baryon acoustic oscillations, etc.) indicate in favor of vacuum 
energy /ACDM cosmology {w = —1). However the ultimate resolution of question about nature of dark energy seems 
to be impossible due to sufficiently large uncertainties in the determination of the DE EoS parameter w — — 1.04j;g''j'o 

The phantom model {w < —1) firstly proposed in ref.[8| probably is the most exotic model of dark energy. This 
model does not contradict the cosmological tests based on present data although from theoretical point of view the 
phantom field is unstable @ because the violation of all energy conditions occurs. 

The simplest phantom energy model with w = const leads to so-called Big Rip future singularity d, [Tol - [l3j . The 
scale factor becomes infinite at a finite time. Another types of finite-time future singularities are type II singularity [l^ 
where the second derivative of scale factor becomes infinite at finite time while the first derivative is finite and type 
III singularity where the first derivative of scale factor diverges. 

If w decreases sufficiently rapidly with increase of phantom energy density the so-called Little Rip or Pseudo-Rip 
may be realized (l"5l - [T7| . In the first case the acceleration of universe increases but the scale factor remains finite 
always. Of course, the disintegration of bound structures occurs as in a case of Big Rip cosmology. For Pseudo-Rip 
universe, dark fluid energy-density asymptotically tends to constant value, i.e. this phantom energy model mimics 
vacuum energy from some moment of time. 
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To describe dark energy fluid we use the equation of state (EoS): 

PD = F{pd) , (2) 

where F is a function of the energy density. The evolution of the universe then depends on the choice of the EoS. 

In this paper we confront DE models mimicking currently ACDM cosmology with the combined observational 
data including the luminosity distance modulus vs redshift for SNe la, the baryon acoustic oscillations (BAO) and 
matter and DE density perturbations. We demonstrate that the models under discussion are totally viable and not- 
distinguishable from ACDM, while their future evolution may vary in a number of ways. The paper is organized as 
follows. In section II we present brief overview of the EoS fluid formalism. In section III the main constraints from 
observational data are discussed. In section IV we confront DE models with Big Rip or type III future singularity 
as well as Little Rip model with SNe data, BAO data and matter perturbations. It is established the excellent 
coincidence with ACDM model predictions at current universe. The consideration of DE and matter perturbations 
shows that growth of the energy-density may occur at sufficiently small background density but still before the possible 
disintegration of bound objects in the Rip universe. This new effect may bring the future universe to the separation 
in the several domains with different values of energy-density, i.e. kind of chaotic universe may emerge. Section V 
is devoted to comparison of DE models of Pseudo-Rip type or with type II future singularity with combined data. 
The region of parameters where such models effectively coincide with current ACDM cosmology are established. In 
Conclusion section some outlook is given. In the Appendix we briefly describe the properties of Quasi-Rip universe. 



II. SCALAR DARK ENERGY MODELS 



For the spatially-flat Friedmann-Robcrtson- Walker (FRW) universe with metric 

ds^ = dt^ -a^{t)(dx^ + dt/ + dz^), (3) 



the FRW equations are given by 



where p and p are the total energy density and pressure, a is the scale factor, ' — d/dt, and we use the natural system 
of units in which SttG = c = 1. For simplicity we assume that universe is filled only by cold dark matter, baryon 
matter and dark energy i.e. p = po + Pm and p ~ pij. 

One can rewrite the dark energy EoS ([2|) in the following form: 

PD = -pD -5(pd), (5) 

where g^pu) is some function. The case (?(pd) > corresponds to the EoS parameter w < — 1 (phantom) while the 
case gipo) < corresponds to the EoS parameter w > —1. Assuming that dark energy dominates, one can neglect 
the contribution of other components (matter, dark matter). Then from Eq. (jj]), one can get the following expression 
for time variable: 

1 r° dpD 
V o J POO pj^ 9[Pd) 

For current time we choose tg = 0. The quintessence energy-density decreases with time {pu < Pdo)^ while the 
phantom energy-density increases (pu > poo)- Fot' scale factor we have the following expression: 

1 dpD \ 



For simplicity oq = 1 is chosen. A finite-time singularity occurs if the integral ([5]) converges at p — >■ oo. The scale 
factor in this case may become infinite (Big Rip) or may remain finite but a singularity (p — >■ oo) occurs. This is a 
Type III singularity [l8|. The type II singularity realizes if g(p£)) — >■ ztoo at pD = Pf- The pressure of the dark energy 
becomes infinite at a finite energy density. The second derivative of the scale factor diverges while the first derivative 
remains finite. 

Non-singular evolution corresponds to the case when integral (O diverges at — > oo (Little Rip [l^, [TtI. [22j) 
or ^ diverges at po ^ P/[l5j. In the last case the dark energy density asymptotically tends to a constant value 
("effective cosmological constant" or pseudo-rip). 

The so-called quasi-rip [l9| . poj corresponds to (inverse) crossing of phantom divide line ("(de)phantomization" 
[2l|). One can show that in terms of EoS this case corresponds to multiply-valued function g{pD) and phase transition 
at some value of scale factor may occur. 
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III. OBSERVATIONAL DATA 



Confrontation of the theoretical models with observational data includes the comparison with several observational 
constraints: 

(i) the luminosity distance moduli to type la supcrnovac from the Supernova Cosmology Project (23j . 

(ii) BAO (see for example (U), 

(iii) the data for the growth factor for density perturbations (see, for instance.refs. l25l . (26j ) from Lyman-a forest 
in the Sloan Digital Sky Survey and galaxy redshift distortions (for instance.rcfs.[27l|. |28l|). 

Let us consider these observational bounds in detail. 

(i) The distance modulus for a supernova with redshift z = ao/a — 1 is 

/Lt(2:) = const + 5 log 13(2) , (8) 
where D{z) is the luminosity distance. As is well-known 

Dl{z)^-^{1 + z) f h~\z)dz, h^iz)=piz)/po. (9) 
^0 Jo 

Here c is speed of light and Hq is Hubble parameter. We used Hq = 72 km/s/Mpc from the Hubble Space Telescope 
key project [2§|. The best fit for SNe is given in the framework of ACDM cosmology. For such a model (which we 
call the "standard cosmology" (SC)), one obtains 

h{z) = {fl„M{l + zf + QaoY^^ (10) 

Here, flmO is the fraction of the total density contributed by matter at present time, and VIa is the fraction contributed 
by the vacuum energy density. Following the approach of ref.HO] one can exclude the SNe samples satisfying the 
condition j^ofes — MSG I /fobs > 1-9. Therefore in the data set of 580 SNe from [3l| we exclude the following 50 SNe 
samples: 

1992bs, 1992bp, 1995ac, 1999bm, 1997o, 2001hu, 1998ba, 04Pat, 05Red, 03D4au, 
04D3gt, 04D3cp, 03D4at, OSDlfc, 04D3co, 03D4dy, 04D3oe, 04Dlak, 03Dlco, bOlO, 
1995aq, f076, g050, k430, ml38, 2006br, 2006cm, 2006cf, 2007ca, 2004gc, 
10106, 2005U, 20051p, 2005fp, 2005gr, 2005ia, 1997aj, 1308, el40, d084, 
2002hu, 2002iu, 2003ch, 2005gs, 2005hv, 2005ig, 2005jj, 1997k, gl20, 05Str. 

(ii) In Ref. [2j| it was suggested that instead of taking the position of an acoustic peak one should measure the large- 
scale correlation function at 100/i~^ Mpc separation using the sample of 46748 luminous red galaxies (LRG)selected 
from the SDSS (Sloan Digital Sky Survey) sample. The appropriate quantity to be measured is known as A parameter 
and reads as 

A^Dvizor-^^^ (11) 

ZqC 



In Eq. dTTT) the dilaton scale Dv{z) is defined as 

Dv{z) 



DMizfcz^'"" 



H{z) 



where Dm{z) is the comoving angular diameter distance. The redshift zq = 0.35. We have for the parameter A the 
following relation 



1 i-zo 

— / h-^{z)dz 
zo Jo 



2/3 



(12) 



The observational value of BAO parameter is A = 0.469 ± 0.017. 

(iii) As it was shown in ref. [l^ one can neglect density perturbations of dark energy. In this case the dark matter 
perturbations effectively decouple from DE perturbations. The equation that determines the evolution of the density 
contrast ^ in a flat background filled by the matter with density pm is 

Srn + 2H5m = ^PrnSrn (13) 
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z 


fobs 


Ref. 


0.15 


0.51 ±0.11 


[32j, [33J 


0.35 


0.70 ±0.18 


m 


0.55 


0.75 ±0.18 


[35j 


1.4 


0.90 ± 0.24 


[36] 


3.0 


1.46 ±0.29 


[37] 



Table I: The available data growth function / at various redshifts from the change of power spectrum Ly-a forest data in SDSS. 



It is convenient to introduce the function of growth rate of perturbations / = dlnSm/dlna. Using FRW equations 
one can get the following equation for /: 

'^^ ±,f ±f4 + 2V/-|a„ = 0, (14) 



din a •' \ iJ2 



where 51„j is matter fraction of the total energy-density: 17„j ~ O„io(l ± z) /h (z). Finally, using relation 

(1 + ^)37 



d In a dz 

and taking into account that 



-PD 



(1 + .)f + + + ^f^. + / - = 0. (15) 



we get 



where fljj — po/p- For dark fluid with given EoS one can find DE density as function of redshift z. Then, Eq. (|15p 
can be solved numerically. The observational data for growth function / at various redshifts are given in table I. 
In the data analysis we use statistics. The value for some physical quantity x is given by equation 

Xl = (^^^ ~f (16) 

where Xth is theoretically predicted value of x, Xth is experimentally measured value and ax is standard deviation. 
For the data set the total is the sum of all Xx ■ 



IV. DARK ENERGY MODEL I: LITTLE-RIP, BIG-RIP AND TYPE III FUTURE SINGULARITY 



Let us start from the simple model with following EoS: 



fpn_\ 



giPo) = a Pno (17) 

\PDaJ 

where a and /3 are dimensionless constants. If /3 = 1 we have ordinary phantom energy model with constant EoS 
parameter w = ~1 ~ . From Eqs. Q one can see that for various /3 the model PT)) describes three types of 
future universe evolution: 

(a) Little Rip if /3 < 1/2, 

(b) Big Rip if 1/2 < /? < 1 and 

(c) type III singularity if > 1. 

The time left before finite-time future singularity can be estimated as following 
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DE models with type III singularity: analysis of SNe data 
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DE models with Big Rip singularity: analysis of SNe data 
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0.59 


295.09 


352.03 


0.59 


59.0 


352.01 


0.59 


32.8 


352.09 


0.5 


0.56 


242.31 


354.06 


0.56 


48.0 


354.48 


0.57 


26.7 


354.49 



Table II: The optimal value of the parameter Qdo, time left before future singularity and corresponding value of for some 
and /3. 



Note that poo = ^doPq = SflooHQ and therefore 

2 1 

3a2(2/3- l)r!^'^o 

Let us restrict ourselves to cases (b) and (c). The results for Little Rip models are similar. From Eq. one can 
derive the dependence of dark energy from redshift: 

^^(,) = |pz.o{l-3a2(l-/?)ln(l + z)}^, /3^1, (19) 

[ PD0(l + 2:)-3" , 

SNe data analysis only. Using Eq. ([19]) one can calculate the theoretical dependence fi{z) from ([9]). Of course 
the best-fit of SNc data is ACDM-model {a = 0). For = 0.72 the value of Xsn is minimal: Xsn === 347.06. For 
a and given f3 one can select the parameter floo to describe the experimental data with good accuracy. In table 
II the optimal floo for some values of parameters /3 and a are given. We also calculated the Xsn fo'^ comparison with 
ACDM-model and the time left for future singularity tf (we restrict ourselves to case (3 > 1/2). 

One can see that in principle the parameters of our model can vary in sufficiently wide limits. However, the 
consideration of data for matter perturbations allows to restrict the range of these parameters. 

Analysis of matter density perturbations data. Again the best-fit for observational data is ACDM-model. 
For Qdo = 0.72 we find that Xdp + Xsn is minimal at /(O) = 0.503. The value of Xdp is 0.879 (within 95% C.L). The 
consideration of matter perturbations data only leads to minimal value xi)p = 0-34 at ^Iuq = 0.78 and /(O) = 0.441. 
However, in this case the Xsn is sufficiently large: Xsn — 370.01 although this lies within 95% confidence level. 

For given (3 one can further vary parameter a and find the optimal values of ^do for fitting SNe data. Then one 
can find /(O) from approximation of matter perturbations data. Hence, we find the maximal value of a at which our 
model fits the matter perturbations data with 95% C.L. The results of our calculation for the model with future Big 
Rip singularity are given in table HI. For DE model with type HI future singularity the results are similar (see also 
[40|). We find that maximal values of parameters a'^^^, flno and /(O) in fact do not depend from 13. It is interesting 
to note that SNe data description for these models coincides with ACDM cosmology up to excellent accuracy. 

The current EoS parameter for the model under consideration is 

Wo = —1 — . 

Therefore small values of parameter correspond to small deviation of w from —1. 

The dependence of In a as function of the time is depicted for some parameters (figure 1). For the illustration, on 
the next figure the comparison between observational data for density perturbations and theoretical predictions is 
made. There is no significant difference between standard cosmology and our models. 
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DE models with Big Rip or type II future singularity: analysis of SNe+DP data 





0.55 


0.75 


0.95 


1.5 


2 


3 


2 

a. 

^max 


0.03 




0.71 


/(O) 


0.514 


tfmin, Gyr 


3.6 X 10^ 


717.3 


398.5 


179.3 


119.6 


71.7 


XSN 


347.28 


A 
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Table III: The maximal value of parameter (and correspondingly the minimal time before singularity) for given /3 at which 
the model (|17p describes the density perturbations data with 95 % C.L. The optimal value of Qdo for fitting SNe data, Xsiv 
and BAO parameter are also given. 



14 
12 



10 



Ina 8 



6 



4 



2 





20 40 60 80 100 120 140 160 180 200 

t, Gyr 

Figure 1: The dependence of scale factor from the time for model (|17p with = 0.02, /3 = 0.95 (big rip) and 13 — 3 (type III 
singularity). In the interval ~ 50 Gyr the universe expansion in these models is very close to the one expected from ACDM 
cosmology. 

Another DE model with type III singularity. One can consider another DE model which is more close to 
ACDM-cosmology. This model firstly was considered in [s^, [s^. From the following EoS 

F{p) = -(3'a}p'j^'/' , (20) 

where af, and e are positive constants one can find the dependence of the dark energy density on the scale factor 

PD = r'^^Hal - a')-^/' . (21) 
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Figure 2: The comparison of observational data with theoretical predictions for the model (jlTp with o? = 0.02, /3 = 0.95 (Big 
Rip) and /3 = 3 (type III singularity). One can see the coincidence with ACDM cosmology. 



For dimensionlcss Hubble parameter as function of rcdshift we have therefore 



h^{z) = f}„(l + zf + fiD(l + zf 



iVo(l 



- 1 



3/e 



(22) 



One can see that for large iVo our model mimics ACDM cosmology with excellent precision. Therefore our model 
can fit the Supernova Cosmological Project data. For iVg ^ 1 the dark energy density is nearly constant in the 
interval < t < i-e. the model (PU)) mimics a cosmological constant in the past but it leads to a finite-time future 
singularity. 

The current EoS parameter is 



Wo 



iVo- 1 



For given value of and e one can find that such model describes the observational data with good accuracy. The 
results are given in table IV. One can see that model ([20]) is more close to ACDM cosmology than DE model with type 
III singularity considered above. The model describes supernova data with the same precision as the ACDM model 
although for density perturbations data the agreement is slightly worse than in the case of standard cosmological 
model. It is interesting to note that agreement between this DE model and observational data is better for large 
values of e. The model (|20p coincides with ACDM model in past with excellent precision but its future evolution 
shows radically different dynamics (see fig. 3). 

The future evolution of dark energy and matter density perturbations. To conclude this section we 
consider the question about evolution of matter density perturbations in future. In such a case one should account 
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tf, Gyr 


13.5 
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3.2 


18.1 
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4.1 


24.7 


10.7 


5.4 



Table IV: The optimal parameters Q,do and /(O) for fitting SNe and density perturbations data for model (|20|l for various wo 
and e. The time before future singularity is given also. The coincidence with ACDM model is better for larger values of e and 
smaller values of wo. 
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Figure 3: The dependence of scale factor from time for DE model (|20|l with various parameters choice. For comparison the 
same dependence for ACDM cosmology (straight line) is given. 



the dark energy density perturbations. The complete system of cosmological perturbations equations can be found 
for example in refs.ji^l, Ell- the system of units in which SttG = c = 1 and for the case when effective sound speed 
for dark energy is equal to 1 these equations can be written in the following form 

S'lj + 3(1 - WD)aH6D - (1 + wd)S„. + (1 + wd) (k + 9a^H'^^-—^) Vd = 0, (23) 



- 2aHVD - 5d = 0, (24) 

l + WD 
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C + —S'„i - liPmSrn + (1 + Swd)pd5d) = 0. (25) 
a 2 

Here denotes derivative on conformal time t {t ^ J a^^dt), SD = Spo/ PDi Vd is the velocity perturbation of dark 
energy and = dpn/dpo is the adiabatic speed of sound. The fc is a wavenumber of the corresponding mode. One 
can rewrite Eqs. ([23 |) -(f25 |) using the relation 

da 

dr = 



a^H 



Thus, we have 



^ + -(1 - wd)Sd - (1 + wd)^ + 4r(l + ^d) (~k + 9a'h'^-S ]Vd = 0, (26) 
da a da a'^ti \ k 

l^,„.„, ,27, 

da a n 1 + wd 

d'^Srn . /3 . dh\ dS,n 1 1 / r , n 1 Q ^ X \ t^Q\ 

(PmOm + (1 + 3w_D)p£)dD) , (28) 



da'^ V a hda J da 2a^h'^ H, 



where h = H/Hq is dimensionless Hubble parameter and k is dimensionless wave number measured in the units of 
Hq/c. One notes that last equation of this system is equivalent to (jl3p if the dark energy density perturbations are 
neglected. However, the future DE density increases and even its small perturbations may become important for the 
analysis of matter density perturbations (see r.h.s. of Eq. 

For the model ([TT]) the EoS parameter and adiabatic speed of sound arc given as functions of scale factor 

w{a) ^-l~a^ (1 + 3a2(l -;3)lna)"^ , = -1 + /?(! + w(a)). 

For current value of the scale factor we put simply ag = 1. Initial conditions for integration are given for example 
in ref.^2|- The main results of numerical analysis of system ((26)) - ((28|) arc the following. For DE model p7|) which 
evolves to Big Rip at acceptable values of parameters (see Table 3) at small scales (large wavenumbers) the amplitudes 
of matter density perturbations increase insignificantly in comparison with 5mo and remain constant. The same 
picture takes place in ACDM cosmology. However, for the standard cosmological model all modes of matter density 
perturbations evolve in the same way. In the case under consideration for sufficiently large scales (> 2000 Mpc) the 
decay of matter perturbations occurs. This decay occurs faster for perturbations with larger scales (see Fig. 4). The 
value of dark energy perturbations grows very quickly before a ~ 10 but then increases very slowly. If /? > 1 (type HI 
singularity) the picture looks the same for matter density perturbations but character of dark energy perturbations 
evolution changes. We observe the approximately linear growth of 5d as function of In a in the interval 1 < a < 10 
but then 5m increases faster according to the law Sq ^ (ln(a//a))~^. Shortly before the singularity the perturbations 
become large. For example for /? = 3, a'^ = 0.03 the moment of singularity corresponds to a/ = 258.67 {tf = 71.7 Gyr) 
and (5£) ~ 1 for fc = 0.0005 Mpc~^ at a ~ 235 [t = 70.4 Gyr). It is interesting to note that growth of density contrast 
occurs at sufficiently small background density pD < lOpjyQ. The sharp growth of density perturbations occurs before 
possible disintegration of such bound structures as Solar System or Milky Way due to enormous acceleration of 
universe. One can assume that like the formation of large-scale structure due to matter density perturbations in the 
early universe this growth can lead to the formation of separate domains with various values of dark energy density 
( "dark large-scale structure" ) . The cosmological dynamics within these " sub-universes" is defined by the background 
dark energy density. In fact, we cannot speak about the uniform evolution of such universe. It is possible that in some 
domains the perturbations lead to the decrease of DE density and induce the de-phantomization (one can assume for 
example that EoS in the form pT| is valid only for pn > p*j~, but if p < the phase transition occurs such that 
w > —1). Yet, in another domains the perturbations of DE energy-density lead to singularity faster than it would be 
expected without consideration of density perturbations evolution. In a sense this picture maybe similar to chaotic 
inflationj4l| reversed in time. The homogeneous universe ends its existence in "chaotic" state. Of course, this is 
rather speculative possibility. 

V. DE MODEL II: PSEUDO-RIP AND TYPE II SINGULARITY 

Let us consider another DE model with EoS 

g{p)^a^PDo[l-^y , 0<PD<Pf, f3^0. (29) 
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Figure 4: The evolution of matter density perturbations (top) and DE energy-density perturbations (bottom) in the model 
(|17|l for the cases with Big Rip (bottom lines) and type III singularity (dotted lines). For DE perturbations the evolution of 
the mode with k = 0.0005 Mpc~^ is depicted. For matter density perturbations the evolution of modes with k = 0.05 Mpc~^ 
(1) and k = 0.0005 Mpc~^ (2) is shown. The model parameters are chosen as = 0.03 and /3 = 0.95 and /3 = 3. 

Such dark fluid leads to the following variants of evolution: 

(i) DE energy-density asymptotically tends to p/ if /3 > 1/2. Therefore, the universe expands according to de Sitter 
law at t oo (Pseudo-Rip). 

(ii) DE energy-density reaches pj for tj <oo ifO</3< 1/2. 

(iii) type II singularity occurs if /3 < 0. The second derivative of scale factor diverges while first derivative remains 
finite. 

Let us consider two simplest cases: (3 = 1, —1. 

(a) /3 = 1. For EoS (|29|) the dark energy density as function of redshift is 

Pd{z) = PuoA-i (l - (1 - A)(l + z)3"'^) , A - pDo/Pf. (30) 

As in previous section we consider only SNe data first. For various A one can find optimal values of parameters a 
and ri_Do- Results are given in table V. 

The analysis of density perturbations data shows that only in the narrow range of these data can be described 
with 95% C.L. (see table VI). This fact simply means that current value of EoS parameter 

Wo = -1- a^{l- A) 
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DE models with Pseudo-Rip: analysis of SNe data 





A 


= 0.1 


A 


= 0.2 


A 


= 0.5 


A 


= 0.8 


A = 


= 0.95 






Xsiv 




XSN 




XSN 




Xsiv 




Xliv 


0.1 


0.68 


347.30 


0.68 


347.25 


0.69 


347.15 


0.71 


347.04 


0.71 


347.00 


0.2 


0.64 


348.17 


0.65 


347.85 


0.67 


347.38 


0.70 


347.07 


0.71 


346.96 


0.3 


0.61 


349.17 


0.62 


348.72 


0.65 


347.71 


0.69 


347.12 


0.71 


346.99 


0.4 


0.59 


350.54 


0.60 


349.83 


0.63 


348.10 


0.68 


347.20 


0.71 


347.09 


0.5 


0.56 


351.66 


0.57 


350.69 


0.61 


348.50 


0.67 


347.28 


0.70 


347.00 



Table V; The optimal value of parameter and corresponding value of Xsjv for some o? and A. 
DE models with Pseudo-Rip: analysis of SNe-|-DP data 



A 


0.1 


0.2 


0.5 


0.8 


0.95 


^max 


0.03 


0.03 


0.06 


0.13 


0.45 




0.71 


/(O) 


0.514 


Xsiv 


347.19 


347.11 


347.30 


347.19 


346.98 


A 


0.482 



Table VI: The maximal value of o? for various A and /(O) at which the model p9|l describes the density perturbations data 
with 95% C.L. The optimal value of Q,dq for fitting SNe data, xliv and BAO parameter are also given. 



only insignificantly deviates from —1. 

Note that same analysis can be performed for quintessence DE model with asymptotic de Sitter evolution. The 
parameter A in this case is larger than 1. The Eq. pop remains correct. The analysis shows that for A '--^ 1 and small 
o? the model describes observational data with good precision. 

(b) /3 = — 1 . In this case for DE energy-density as function of rcdshift wc can derive the relation 

Pd[,z) = pz5oA-i {l - ((1 - A)2 + 6a2Aln(z + l))i/2| (31) 
The time left for future singularity can be found from ([5]) by integrating from p^io to pf. We have 

Next, as in the previous cases we calculated the optimal value of for fitting SNe fata for given A and (table 
VII). One sees that only at sufficiently large o? and A our model declines from observational data significantly. For 
A 1 the model (|29p fits the observational data only at — > 0. 

For current value of EoS parameter in the model with type II singularity we have 



1- A' 



Again, wq is close to —1. 

The combined analysis of SNe data and data for density perturbations gives the results similar to previous cases. 
The deviation from ACDM model is sufficiently small for fitting observational data although these data in principle 
do not prevent the possibility of future singularity within ^ 100 Gyr (table VIII). 

The analysis of the evolution of density perturbations can be performed as at the end of previous section. Let 
us briefly recall the main results. For Pseudo-Rip the matter density perturbations evolve in the same way as in 
ACDM model. For dark energy density perturbations we have rapid growth before a '-^ 10 and then slow decay. For 
models with type II singularity the picture coincides with that of ([TT]) with type III singularity. Hence, DE models 
mimicking ACDM and fitting current observational bounds may show different exotic behaviour in the future: finite- 
time singularities, disintegration of bound structures (Little Rip or Pseudo-Rip cosmologies) or decay of cosmological 
perturbations. 
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DE models with type II singularity: analysis of SNe data 





A = 0.2 


A = 0.5 






i/, Gyr 


Xliv 




i/, Gyr 


xIat 


0.1 


0.67 


27.6 


347.67 


0.64 


8.7 


348.64 


0.2 


0.62 


14.2 


349.02 


0.59 


4.6 


351.75 


0.3 


0.59 


9.8 


351.04 


0.54 


3.2 


356.14 


0.4 


0.55 


7.6 


353.18 


0.51 


2.5 


360.67 


0.5 


0.52 


6.2 


355.65 


0.48 


2.0 


365.51 



Table VII: The optimal value of parameter ilDo, time before singularity and corresponding value of Xsiv for some c? and A. 
DE models with type II singularity: analysis of SNe+DP data 



A 


0.2 


0.5 


^max 


0.02 


0.012 




0.71 


tfmin, Gyr 


134.1 


69.4 


/(O) 


0.513 


2 

XSN 


347.13 


347.10 


A 


0.482 



Table VIII: The maximal value of for various A and /(O) at which the model 1)29^ describes the density perturbations data 
with 2a CL. The optimal value of Qdo for fitting SNe data, Xsn a-nd BAO parameter are also given. 



VI. CONCLUSION 



In summary, wc confronted number of DE models mimicking ACDM epoch with current EoS parameter being 
very close to —1 with combined observational data: SNe data, baryon acoustic oscillations data and DE and matter 
energy-density perturbations. It is explicitly demonstrated that there exists sufficiently wide region of parameters for 
each of DE model under discussion where these theories are not less viable than the standard ACDM model. On the 
same time, DE models under consideration show qualitatively different future behaviour with Big Rip, type II and 
type III future singularities. Little Rip, Pseudo-Rip or Quasi-Rip evolution. Nevertheless, current observational data 
cannot determine whether or not the universe will end in a future singularity. 

It should be noted that the observational data for density perturbations are more sensitive to the deviation from 
standard cosmological model in comparison with SNe data. Moreover, the account of density perturbations in Rip 
cosmology indicates to sharp growth of density at sufficiently small background density still before the possible 
disintegration of bound objects. This growth from our viewpoint can lead to possibility that future universe may 
split in the number of separate regions so that it becomes chaotic and never reaches the Rip singularity. Further 
consequences of above effect will be considered elsewhere. 
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Appendix: DE model with dephantomization: Quasi-Rip 

In this Appendix we briefly discuss the possibility of so-called "Quasi-Rip" scenario considered in ref.[l^. Let us 
consider the case when DE pressure depends from its energy-density as follows: 



-PD 



PD 



2 9 



-pD + ^a^PD 



hlPlIL 
PD 



1/2 
1/2 



a < ax, 
a > ar, 



(33) 
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Figure 5: The pressure as function of energy-density for model (|33p . The dotted line corresponds to vacuum energy with 
w = —1. 

where po < Pm and Inax ~ ot'^ (lii(pm/pz3o))^^^ is a value of scale factor at which the dephantomization occurs. In 
this moment the value of DE energy-density reaches the maximal value pm . For small o? the universe acceleration is 
maximal for pD k, p,„ exp(— a^/2). The pressure as function of energy-density is depicted on Fig. 5. 
The dependence of DE energy-density from redshift z for EoS ([55]) is 



The Eq. ([55]) coincides with the equation considered in ref. [I^. As is demonstrated in above work this model fits 
modern SNe data in the large interval of free parameters. One can show that the model p3p is consistent with matter 
density perturbations data also. Hence, the Quasi-Rip model maybe consistent with current observational data. 

The specific feature of Eq. ([33)l is that scale factor is double- valued function of the energy-density. The branch point 
corresponds to dephantomization: the pressure in this point is defined uniquely. It is obvious that the described scheme 
can be generalized. For equations of state with branch points the evolution of universe contains (de)phantomization 
and Quasi-Rip epochs. 
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